A stringent constraint on neutrino Lorentz invariance violation 
from the two IceCube PeV neutrinos 



Enrico Borriello, Sovan Chakraborty, and Alessandro Mirizzi 
77" Institut fur Theoretische Physik, Universitat Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 



LAPTh, Univ. 



o 

(N 



m 

(N 

X 

43 

Q-i! 

6 



> 
m 

oo 
in 

cn 
o 



is 



Pasquale Dario Serpico 
Savoie, CNRS, B.P.110, Annecy-le-Vieux F-74941, France 



It has been speculated that Lorentz invariance violation (LIV) might be generated by quantum- 
gravity (QG) effects. As a consequence, particles may not travel at the universal speed of light. In 
particular, superluminal extragalactic neutrinos would rapidly lose energy via the bremssthralung of 
electron-positron pairs (v —¥ v e + e _ ), damping their initial energy into electromagnetic cascades, a 
figure constrained by Fermi-LAT data. We show that the two cascade neutrino events with energies 
around I PeV recently detected by IceCube — if attributed to extragalactic diffuse events, as it 
appears likely — can place the strongest bound on LIV in the neutrino sector, namely S = (v 2 — 1) < 
e>(lCr ls ), corresponding to a QG scale M QG > 10 5 M P i (M QG > 1CT 4 Mpi) for a linear (quadratic) 
LIV, at least for models inducing superluminal neutrino effects (S > 0). 



PACS numbers: 11.30.Cp, 95.85.Ry 



I. INTRODUCTION 

It is conceivable that Lorentz invariance might be vi- 
olated in a candidate theory of Quantum Gravity (QG), 
see [l[ for a review. In this context, the space-time foam 
due to QG fluctuations might cause e.g. highly boosted 
energetic particles to propagate at speed v different from 
the velocity of light c. The resulting Lorentz invari- 
ance violation (LIV) effect can be phenomenologically 
parametrized in terms of 5 defined as (we set c = 1) 



dE 
dp 



V 



(i) 



where, assuming that there is at least one frame in which 
space and time translations and spatial rotations are ex- 
act symmetries (typically the lab one), one writes 



E 2 



p 2 + m 2 + f(p, ...). 



(2) 



The function / is often expressed in powers of momentum 
over the phenomenological scale Mqg when LIV becomes 
important, of the type oc p e Mqq . Hence one can trans- 
late a constraint on 6 into a constraint on Mqg via a 
relation of the kind 



vodf_ 
E dp 



\M QG J 



(3) 



with e.g. cubic (£ = 3) or quartic (£ = 4) terms in powers 
of p added to the dispersion relation of Eq. (|3]), inducing 
"linear" (n = 1) or "quadratic" (n = 2) deviations from 
LI occurring at a mass scale Mqg- The + (— ) implies 
superluminal (subluminal) velocity. 

The observation of LIV effects of this nature, if sup- 
pressed by the Planck scale Mpi = 1.22 x 10 28 eV, is 
not possible in any Earth-based experiments, nor it is 
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anticipated that the needed sensitivity will make them 
accessible in the near future. Presently only high-energy 
astrophysics experiments can probe the consequences of 
tiny LIV 0, [|[ . Tests of LIV in photon propagation from 
objects like Gamma Ray Bursts (GRB) |3|, Active Galac- 
tic Nuclei (AGN) d and Pulsars [j| looking into delays 
in arrival times of photons of different energies can probe 
a quantum gravity mass scale up to Mqg ~ 10 26 eV for 
the linear case, and Mqg ~ 10 19 eV for the quadratic 
case. 

Experimental probes of LIV for neutrinos are limited 
by the scarcity of neutrino data from distant astrophysi- 
cal sources. In particular, from SN1987A data, exploiting 
the flight delay of the 7 with respect to v of a few hours, it 
has been obtained a bound of 8 < 10~ 12 Jf|, correspond- 
ing to Mqg > 10 18 eV for n = 1, or Mqg > 10 12 eV 
for n = 2, taking E ~ 10 MeV as typical SN v energy. 
The supernova limit on LIV can reach Mqg ~ 10 21 eV 
(2 x 10 14 eV) for the linear (quadratic) energy depen- 
dence, observing narrow time structure in the neutrino 
emission, like the ones associated with the neutroniza- 
tion burst 0- The bounds on LIV could be further 
improved in case of detection of high-energy neutrinos 
(E > O(GeV)) from astrophysical sources, like Gamma 
Ray Bursts or Active Galactic Nuclei 0]. Indeed, from 
the energy dependent time of flight delay of neutrinos and 
the corresponding 7 rays, induced by LIV, the limit that 
would be placed is Mqq ^ 10 26 eV for the n = 1 case, 
and Mqg > 10 19 eV for the n = 2 case, respectively. 

For the specific case of superluminal neutrinos, LIV 
would also allow processes that would be otherwise 
kincmatically forbidden in vacuum, namely neutrino 
Cherenkov radiation [y — > ^7), neutrino splitting [y — > 
vvv) and bremssthralung of electron-positron pairs 
(y — > v e + e~) (see, e.g. All these processes would 
produce a depletion of the high-energy neutrino fluxes 
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during their propagation. In particular, among the pre- 
vious processes, the neutrino pair production has been 
recognized as the fastest energy-loss process for LIV neu- 
trinos and has been used in jiol [Tlj to invalidate the su- 
perluminal velocity claim made by the OPERA Collabo- 
ration [l2^ . A similar analysis has also been able to put 
strong constraints on supcrluminal velocities of higher 
energy neutrinos. The observation of upward going at- 
mospheric neutrino showers (with path-length L ~ 500 
km), measured at E > 100 TeV at the IceCube exper- 
iment, has allowed to put the strongest bound on LIV, 
namely S < 10 -13 [l3[. In Sec. II we explore a new strin- 
gent bound based on this mechanism, following from re- 
cent observations of IceCube. In Sec. Ill we discuss these 
results and conclude. 



II. A NEW STRINGENT BOUND 

Since the effects of LIV processes typically grow with 
energy and distance, it is worthwhile to exploit high- 
energy neutrinos traveling along a sufficiently long base- 
line. It has been argued for example that observations of 
the cosmogcnic neutrino background (produced by ultra- 
high energy cosmic ray losses) might lead to stringent 
constraints, see [lij]. In this context, the IceCube exper- 
iment has recently reported the detection of two cascade 
neutrino events with energies around E ~ 1 PeV in the 
period of data taking 2011-2012 At the moment, the 
origin of these events is not settled. It seems that conven- 
tional atmospheric neutrinos or cosmogcnic neutrinos [l6[ 
arc unlikely to be the sources of these events, and that 
the chance that these events are associated to prompt at- 
mospheric neutrinos from heavy quark decays is at most 
a few % [TtJ • There is also another intriguing experimen- 
tal indication: an IceCube analysis of neutrino-induced 
muon track events shows that, although the sample is 
dominated by conventional atmospheric neutrinos, data 
do prefer marginally (at the 1.8 a level) an extra com- 
ponent at E > 100 TeV [Hj|. Albeit of little significance 
at present, it is intriguing that the extrapolation of this 
best-fit flux to PeV scale is in perfect agreement with 
the observed two PeV cascade events mentioned above. 
Although the significance of an additional component is 
still low, the hypothesis that these two neutrinos could be 
the first indication of a diffuse extragalactic astrophysical 
flux at the PeV scale appears at the same time the most 
likely and the most exciting one. 

In this section, we deduce an important improvement 
over the current bounds on supcrluminal neutrino veloc- 
ity, following from the above-mentioned assumption. In 
Sec. IIIII we shall comment further on this point. 

First, let us remind that for S > the process v — > 
v e + e~ is kincmatically allowed provided that 

E v > 2 m e /y/8 ~ PcVVlO" 18 ^ (4) 

(and assuming LI conservation in the electron sec- 
tor) Q3 . This process implies a pair production decay 



rate of 0: 

i ri2 p5<>3 

re± = nill^ = 2 - 55xl053<53 ^ vMpc_1 - (5) 

As in [3, we neglect here the process (y — s- vvv), which 
is only expected to bring minor modification to our ar- 
gument. On the other hand, if the process associated 
with the rate of Eq. ([5]) is forbidden because of threshold 
effects, the v — > wy is anyway operational and a channel 
for energy losses, although it is two to three orders of 
magnitude less efficient than v — > ve + e~ . 

Remarkably, without any significant additional infor- 
mation on astrophysical sources of neutrinos (neither on 
the specific mechanism of neutrino production), this is 
sufficient to derive a strong bound. Let us assume that 
the PeV neutrinos are generated somehow in extragalac- 
tic sources. If the above processes are effective, the neu- 
trino flux will be soon depleted at the expense of injecting 
electromagnetic energy while propagating to the Earth. 
The interaction length of e ± onto CMB photons is ex- 
tremely short, of the order of few kpc. Even a gamma 
ray will not propagate much more than O(10) kpc at 
these energies before pair-producing on the CMB (see 
for example [191]). The inverse- Compton scattering pho- 
tons induced by e^ will accumulate a gamma-ray flux 
at energies below 0(100) GeV. This flux is constrained 
by Fermi data [20| to have an integrated energy den- 
sity u^= (4tt/c) / Edip^/dE < 5.8 x 10 _7 eV/cm 3 , see 
e-g- [2l|. Now, the two events detected in IceCube in 
two years of data taking would imply a diffuse energy 
flux (E^d(p E /dE) ~ 10~ 8 GcV enr 2 s" 1 sr" 1 0, actu- 
ally quite close to the so-called Waxman-Bahcall bench- 
mark [22I ] . These numbers imply that one cannot tolerate 
an energy density wjf in the initial neutrino flux larger by 
about two orders of magnitude large than the observed 
one w° bs , otherwise the electromagnetic energy injected 
via the process would basically saturate the bound. Once 
accounting for the fact that the process in question trans- 
fers a lar ge f raction of the initial neutrino energy into the 
e^ pair [id ], a simple back-of-the envelope calculation 
leads to an approximate constraint of the type: 

, .obs 

e -">^.io-2, (6 ) 

where d is the characteristic distance of the sources. Pro- 
vided that the channel v — > v e + e~ is open, we get 

S 3 e ± d Mpc < 1.8 x 10~ 53 . (7) 

For the simplest scenario of cosmologically distant 
sources, a reasonable value is d ~ O(10 3 ) Mpc. Equa- 
tion (0 nominally implies 8 < O(10 -19 ), which means 
that the pair-production mechanism is not operational 
and that the actual bound is thus 5 < O(10~ 18 ) (see the 
condition of Eq. (0J). A priori, the v — > ^7 is kinemati- 
cally accessible to PeV neutrinos even for 5 < C(10~ 18 ) 
and depending on its rate it could put a more stringent 
bound, but it is easy to check that it is not the case: it 
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leads instead to a slightly weaker (albeit of comparable 
order of magnitude) bound. 



III. DISCUSSION AND CONCLUSION 

We have derived a very stringent bound on LIV in the 
neutrino sector, 5 < 0(1O~ 18 ), from the observations of 
two PeV scale IceCube events and remarkably few as- 
sumptions on the underlying astrophysical sources. Let 
us discuss the model dependence in more detail. We just 
assumed that the events are due to an extragalactic flux, 
hence with characteristic distance of source(s) of the or- 
der of a Gpc. Of course, once additional information will 
be available (e.g. on the number density and redshift 
distribution of the sources contributing to the flux) an 
improved calculation will be possible. We do not expect 
however that the ignorance of these details affects the 
bound very much. In fact, within our approximations 
the bound on S scales as 



S < 1CT 18 x Max 



l.G 



ln(u; 7 /< s ) 

^PoV <^Mpc 



E 



PcV 



, (8) 



with a very weak dependence on the initial assumption 
on the fluxes, as well as a dependence only on the cubic 
root of the typical distance. For example, it is worth 
noting that even in the extreme (and unlikely) case where 
the observed neutrino flux were dominated by the closest 
extragalactic source candidates, at distances of the order 
of several Mpc — Centaurus A, which is the closest one, 
is at about 4 Mpc from us — one would still deduce 8 < 
0(1O -18 ). Again, this is basically suggesting that the 
energy loss channel via pair emission must be a closed 
channel. On the other hand, detecting higher energy 
events would improve and make the bounds more robust. 

The extragalactic nature is essential in making the 
electromagnetic "cascade" argument operational, but ac- 
tually it is also the most likely explanation, given that the 
flux level required to match observations appear close to 
predicted diffuse fluxes (sec e.g. 23]). Not only that: an 
indirect argument against a closer astrophysical origin is 
that the Galactic diffuse flux at PeV energy can be com- 
puted quite reliably and is expected to be much smaller, 
see e.g. [24|. This is true provided that the source is rela- 
tively "thin" to cosmic rays, so that the flux is dominated 



by cosmic ray production in interstellar medium. This 
appears almost unavoidable, especially since we would 
need accelerators capable of maximal energies up to 10 16 
eV/nucleon at least, which would be hard to impossi- 
ble to achieve with large energy losses within the accel- 
erator. Anyway, some bound should exist also if one 
or more hypothetical Galactic sources contribute to this 
flux. Lacking promising candidates of this type in the lit- 
erature, we do not consider this further. In general, the 
bound should be then obtained by estimating an input 
flux upper bound by some other sort of multi-messenger 
constraint: Typically, other by-products (or tertiaries of 
their energy degradation) of the neutrino production pro- 
cesses at the source. 

Finally, what if the events should be eventually at- 
tributed to atmospheric background? We briefly men- 
tioned that this appears unlikely at present. Yet, it 
is worth noting that even in that case, the current 
bound (l3T | could be still improved by a more modest 
amount (perhaps one order of magnitude), although a 
more detailed analysis of the type of [HI — properly ac- 
counting for prompt contribution and its uncertainty — 
would be needed to draw a more quantitative conclusion. 

In summary, we have argued that a confirmation of 
the extragalactic astrophysical nature of the PcV events 
detected by IceCube would not only open a new win- 
dow to the high-energy universe, but also allow a signif- 
icant jump in tests of fundamental physics. To obtain 
this result, as more and more customary in astroparticle 
physics, we relied on a multi-messenger strategy (here the 
link with diffuse gamma ray background), which proves 
once again the power of this approach. 
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